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QHAFTER I 
IHTRODUGTIOIJ 


To the vast field of applioablons of a sonic onduct or 
p— n junction t a significant contribution isr of its voltage 
variable «apacitance f in reverse bias condition, Some of the 
major applications of a p-n junction, as a voltage variable 
capacitance, are in parametric amplification, harmonic 
generation, mixing, detection and voltage variable tuning. 

The field of application of a p-n junction, as a 
varactor (variable reactance) can be divided into two areas. 

(1) where change in capacitance of the p-n junction', with 
voltage is used, as in tuning, 

(2) where the rate of change of capacitance is used as in 
parametric amplification. 

For both the areas of applications, hyperabrupt 
junctions exhibit a 'superior characteristic to that of other 
types of p-n junctions i.e. abrupt and linearly graded junc- 
tions, This is because the former is more voltage sensitive 
and hence provides a larger variation of ccapacitance and 
higher rate of change of capacitanoei The higher voltage 
sensitivity of capacitance in hyper-abrupt junctions is 
attributed to the peculiar impurity distribution near the 
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junction . In view of above mentioned grouping of appli - 
cations, the breakdown voltage (Vg)- and C-V index (n) are the 
important parameters of a varactor, as the former dominantly 
determines the range of variable capacitance available and 
latter, the rate of change of capacitance. Both the 
parameters solely depend upon the impurity distribution on 
both the sides of the junction. 

A. hyperabrupt junction is characterized by the 
uniformly doped one side of the junction, while on the 
other side, the impurity concentration decreases with distance 
from the junction. Uniformly doped side is generally heavily 
doped while <in the other side of the junction, commonly 
known as retrograded region, the impurity distribution may 
be satisfying either of the exponential , complementary error 
function, gaussian function or the power law distribution. 

In one-sided hyperabrupt junctions, therefore, the 
breakdown voltage and C— V index will depend upon the 
impurity distribution in the retrograded region. Since the 
technology, today, permits high controllability of impurity 
distribution, breakdown voltage and C-Y index of a hyperabrupt 
junction can be pre— determined, if their relationship with 
impurity distribution is established. 

Several authors have discussed different aspects' of 
Dreakdown in hyperabrupt junctions [1-5]. Uathanson et.al. 
r4] and M. Shinoda [5] studied avalanche breakdown in silicon 
hyperabrupt junctions. Both the authors, however, in the 
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calculation of breakdown voltage^ have assumed the ionization 
rates of holes and electrons [6-7] to be eq_ual» further 
measurements on silicon p-n junctions by a number of 
investigators^ revealed that this assumption is not correct. 
Ihe ionization rates of electrons were found to be higher 
than that of holes^ in silicon p-n junctions [8-15],# The 
dependence of avalanche breakdown voltage on impurity distri- 
bution, in hyperabrupt junctions, therefore, needs to be 
re-evaluated. Fane is known to have discussed zeneir bre^- 
down in silicon hyperabrupt junctions, 

A number of authors have discussed the C-V characte - 
ristics of hyperabrupt junctions [16-19]. A,Shimiju [16] 
and M, Shi nod a [l7], calculated the 0-V characteristics of 
hyperabrupt junctions assuming exponential and complementary 
error function impurity distributions in the retrograded 
region, respectively. The latter author has also discussed 
the temperature dependence of C-Y characteristics of 
hyperabrupt junctions. The C-Y relationships, assuming 
power law distribution in retrograded region has also been 
derived [18], 

C-Y index (n)j- The 0-Y index is defined as the slope of the 
curve In C vs. In Y where 'C’ and ’Y’ are the capacitance of 
aiid voltage applied across the junction. It can be expressed 
as, 

a(iii c) 


n 


a(ln Y) 


( 1 . 1 ) 
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A comparative study of C-7 characteristics of the 

different types of hyperahrupt junctions reveals that 0-7 

index^{n) varies with applied voltage except in the case 

region 

where the retrograded^has a power law distrilbution. Por the 
junctions in which ’n* varies with voltage, ’n’ is found to 
he first increasing, reaching a maximum and then decreasing 
with increasing voltage# These types of the junctions are of 
practical interest because they can be fabricated by most 
widely used diffusion techniques,. A rigorous study of 
dependence of *n’ on the applied voltage and parameters of 
impurity profile in retrograded region, is of practical 
importance. 

The present work was started with the following 
objectives, 

(1) To, find out the dependence of avalanche breakdown 

voltage (Vg) and space charge width at breakdown 

(W-n), on the parameters which characterize the impurity 
ii 

distribution in the retrograded region. Also, to 
find out empirical relationships, expresssing Vg and 
as a function of the impurity distribution 
parget ers. 

(2) To find out the dependence of 0-T index on the , 

impurity distribution in retrograded region and to 

establish the empirical relationships expressing 

n , the maximum value of C—7 index ,as a function of 
max 

impurity profile parameters. 

* Se£. ^ 
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(3) To fabricate silicon hyperabrupt junctions and 

to investigate how far the experimentally measured 

values of 'V-n* and ’n* are in agreement with the 
b max 

calculat ed results • 



GHAPTBR II 


i 

AVALAKGHE BREAKDOWU II SILICON HYPER ABRUPT JUITOTIOIT 


2.1 Electric field and potential distribution in space 
cbarge region: 

Tile assumed impurity distribution in the hyperabrupt 
junction is shown in Fig.(l). One side of the ju'iption 
has been assumed to be very heavily doped, Qn the other 
side i.e. in the retrograded region, the impurity concentra- 
tion has been assumed to be decreasing exponentially with 
distance from the junction. This choice of impurity 
distribution in retrograded region has been made for two 
reasons. Firstly because it simplifies the calculation 
of the avalanche breakdown voltage and secondly, because 
the exponential distribution has been found to be 'a good 
approximation to the tails of Gaussian arid complementary 
error functions, the types of distributions that are 
obtained from most commonly used diffusion techniques 

Assuming that the impurity concentration decreases 
exponentially with distance from the junction, the 
expression for impurity distribution in retrograded region 
can be written as, 

iKz) = (Kq - Ng) + Fg 


( 2 . 1 ) 



7 





where t 

’x* is the distance from the junction 

is the cross over impurity concentration at the 
junction, on the retrograded re g ion side, 

'!’■ is the characteristic length 

‘Wg* is the background impurity concentration. 

Assuming that there are no mobile carriers in the 
space charge region and also that all the impurity atoms in 
the space cha;t'ge region are ionized, the Poisson’s equation 
for space charge in the retrograded region can be written as? 

■ f e2.2) 

Since the junction has been assumed to be very 

heavily doped on one side, the width of the space charge layer 

nearly all the 

is very small dn this side of the junction and^otential 
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drop occurs across the space charge in retrograded region, 
lEhis yields the following boundary condition for eq,. (2^5). 

V = 0, at s = 0 (2,3)a 

Assuming that the voltage drop across the region, 
outside the space charge layer, is negligible and also that 
the electric field due to inhomogeneous impurity distribution 
in retrograded region is negligibly small compared to the 
field in the space charge region, the second boundary 
condition for eq, (2.2) can be written as, 

E = 0 at x = W {2.3)b 

where ’B’ is electric field 

and V* is width of space charge layer in the retrograded 
region 

fflite integration of eq. (2.2) , with the boundary 
condition (2,3)b yields the equation for electric field 
distribution in the space charge region. 


cs 


B(3c) 


e 


B [ c 1 - 


E 


E 




o 

+ — - — ,(x~¥) ] (2.4) 

^o 

Ihe equation for potential distribution in retrograded 
region can be obtained by integrating eq. (2,4) with 
boundary condition (2.3)a, Tie expression for potential 
distribution can be given as, 
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V(x) = [(1 - 


\ -x/l ' -¥/L 

) • Il» ( 9 ■" ¥* © ) 




+ -2- .(x.w- ^ )I 


(2.5) 


From eqs. (2,4) and (2,5) , the expressions for 
maximum electric field ^ ^ .whicli occurs at x=0 and 

total voltage drop across the junction, ?(¥), can be obtained, 
Ehus, ' 


= - f - [(1 - 


E 


'B 




E 


B 


E 


-¥/L*- ■ 
).( e - 1) 


¥ T 

L 


and 


7{W) = 7^ + V„ = -J- [ ( 1 . ^ ) {1 _ (1- I ) 


Ir 

E, 


( 2 . 6 ) 


-¥/L 

e ) 


E 


B 


E, 


Il-T 

2L^- 


(2.7) 


where, 


2.2 


Y = voltage applied across the junction 
Q» . - 

=3 Built in potential of the junction. 

Avalanche breakdown in a hyper abrupt junction: 

Eie threshold condition for avalance breakdown in a 


p -n junction, where avalanche multiplication is initiated by 


holes, is given by [20], 


% 

a^. exp 
-^0 P 


... X 

'-Jo 


i^r 


a^) dx J: dx = 1 


( 2 . 8 ) 


44. 

and the same for an n -p junction where avalanceh multipli- 
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where, 

’Wg’ is the space charge width in retrograded region 
at breakdown 

’ and ’a^’ are the ionization rates of holes and 
electrons, respectively- 

In the eqs. (2-8) and (2,9), the space charge width in the 
heavily doped region has been assumed to be negligible* 

111 the experimental results have shown that electric 
field dependence of ionization rates can be expressed by 
[ 9 - 10 ] 

a(E) = a. exp [-(b/E)“^] (2.10) 

where 'E* is electric field and’a;b,m are constants* One of 
the early reliable estimates of a,b and m was made by Sze 
and (ribbons [20] from the results of lee et.at [8], More 
recent experiments of Yan Overstraeten and De Man [13] on 
silicon p-n junctions and of ¥ood et, at [14] on silicon 
Schottky barriers have resulted in ionization rates different 
from those obtained by lee et.at [8], In spite of this 
difference, it is interesting to note, that the breakdown 
voltages calculated from these different values of ionization 
rates do not differ from each other to any great extent and 
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are in good agreement with the experimentally measured 
values , both in abrupt as well as linearly graded junctions 
[ 15 ].. ¥e have used the values o.f a,b and m given by Sze 
and Gibbons [20], The values of a and b for electrons, in 
silicon^ are 3*80x10 cm'" and 1*75x10 Y/cm* respectively and 
the same for holes are 2*25x10"^ cm”^ and 3,26x10**^ Y/cm,, 
respectively. ' m is unity for both, the holes and electrons. 

It has been found [29 ] that for the same limits of 
integration, the breakdown field for p’^-n junction is slightly 
higher than that for n‘^-p jimction. However, because p'^-n 

VcSLYtictors 

junctions are invariably used in hyperabrupt junotlone and 
also because the difference is not very signf.fi eani: , the 
breakdown voltaige calculations have been done only for 
p'^-n junctions. The eq. (2,8) has been solved numerically 
for Wj, by Simpson’s rule. 'The sequence of calculation is 
as follows. 

*¥g’' is first given an initial value ’¥’ which is less 

t I 

than the expected value of ¥g. This initial value has been 
selected to be equal to the space charge width at breakdown 
of a one sided abrupt junction, with impurity concentration 
on lightly doped side equal to ’H ’. ’If* is now increased 
in small steps, till reaches ' the value at ^diich 
ionization integral becomes unity. 

T-o determine Vg accurately, the number of subintervals, 
into which the range of integration (0-¥g) has been divided, 
is increased in steps till the increment in the number of 













Nb/Nq =0-10 



Nq ( cm . ’ 
Fig.4{a) 
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subintervals ■ ceases to cause a change in V^. It is relevent 

to point out that the ranges and 0-x have been divided 

into equal number of subintervals. 

lo obtain the breakdown voltage (72)and maximum 

electric field at breakdown been substituted 

for *¥' in equations (2*7) and (2,6) respectively* 

These calculations for'Vg. and have been made for 

a number of values of in the range 3x10 cm to 

17 

5x10 cm with ’1’ given different values in the range 
0.3 micron to 4.16 miron and ^ — given the values 0.1,0.05 
and 0,01, respectively. and have been plotted as 

function of ’ nJ,. with ’ L' as parameter, for the differmt 
values of Eg/'U^ [lig. (2),( 3) and (4)]. 

It can be seen that '7g and Vg very strongly depend 
upon [Fig. (2), (3) and (4)]. ' I’ is seen to affect *7g 

and *¥g, only for lower values of and ’ Ii‘ , Also, a 

comparison of Fig. (2) with Fig, (3) and Fig, (4) reveals , that 

lit! 

Eg/E^ affects 7g and ¥g only for values of ’ L* below 1,06 
micron . The nature of these curves can be explained as 
follows. 

t t II • ‘ 

For the sake of convenience the 7g versus E^ and ¥g 

versus *E plots have been divided into two regions, on the 

1 ! 

basis of value of E^, The region I is that part of the 

I I 17 

plots where E^ 2 x 10 cm and the part of the plots. 

It 17 

where E^ 2 x 10 ' cm ^ is depicted as region II* 
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In region I, for values of ' L' greater than 1.06 

micron it c^an he seen that'w^ is ve3y small compared, to *1*. 

In such a case, the impurity concentration at z=:0,is not 

very different from that at x=¥g ard the junction “behaves 

more or less like an abrupt junction, with impurity 

concentration on lightly doped .side equal to’Kj, The 

calculated breakdown voltages in these cases are nearly same 

as that of the abrupt junctions. 

Ebr those cases in the region I, where 1 ^1.06 micron 

and in region II where 1 ^ 1.06 micron , it can be seen 

that the ratio is much larger than that mentioned in the 

last paragraph. The voltage and space charge width 

for such junctions are significantly greater than those 

for abrupt junction with effective impurity concentration 

equal to These are the cases where hyperabrupt nature 

of the junction significantly shows up. 

In region II for 1 ^ 1.06 micron , the effect of *L * 

and on V ' and V ’ shows up prominantly. In these 

cases, the space charge has covered the major part of the 

retrograded region, before breakdown occurs It is seen 

that Yt, are very sensitive to variations in and 

Nr,/H , in this part of region II, 

The maximum electric field(E^.j^) has been plotted 

against ' nJ, for different values of ’ L’ in Eig. (2)a. 

Because 'E v has been found to be more or less constant with 
mb 

variations in ’SrJ'S » it has not been found necessary to 

JD O 
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plot it for remaining values of expected, 

is found to be maximum for = 5 x 10^"^ cm“^ and is equal 
to 8.64x10^ V/cm, Since tunneling current in silicon p-n 
junctions is significant only at electric fields higher than 
10 V/cm,[2l], the mechanism of breaicdown in all the cases 
considered above, can safely be assumed to be the avalanche 
breakdown , 

2^3 Zener breakdoTO in hyper -abrupt junctions: 

As can be seen from Fig, (2)a, the highest value 
of corresponds to low breakdown voltages (Region I, 

L 1,06 micron ) , As explained earlier the hyperabrupt 
junction. In these cases, behaves more or less like an 
abrupt junction. The zener breakdown which occurs at very high 
electric fields ( y, 10^7/cm,), can be expected to occur only 
in the junctions in which y 5x10 ' cm , and ’ L* is also 
larger than 1,06 micron , so that the space charge width 
is low enou^ to provide high electric field required for 
the tunneling of electrons. Such junctions, however, are 
sure to behave like an abrupt junction, For this reason, it 
has not been necessary to discuss the zener breakdown in 
hyperabrupt junctions. 

tiff 

2,4 Empirical relationship between and 

By linear regression method [22], it has been found 
that and Vg can be approximately expressed as a function 
of and Vg, by the following empirical relations. 
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In (Vg) = 6.92 - 0.85 In - G.II 3 In (I) -0.29 In ( ) 

( 0 . 08 ) ( 0 . 014 ) ( 0 , 016 ) (0.043) ° 

-0.05 In (N^^_^) (2.11) 

(0.04) 

'For the above equation = 0. 973»K^ = 0.972 

and, 

ln(¥g) = 4.63 -1.075 In (F^^) -0.28 ln(l) 

(0.13) (0.025) (0.025) 

- 0.008 (Ng^)"^ 

(0.023) C2.12) 

For equation (2.12), = 0.941* = 0.939 

where. 



— 

R /lO^^ 

on 


0 

^Bn 

= 

R Rg/lO^^ 

1 and 


are in microns 

7g is 

in 

volts 


Equations (2,11) and (2*12) do not hold good in 
the region outside the region hounded hy the chain lines in 
KLg*(2)a and Fig.(2)h, Also the last two terms of the equation 
(2.11) are zero for the cases in which 10 cm Hie 
quantities written below the co-efficients, in brackets, are 
standard errors of respective co-effecients. For difinition 
of R^ and R^ see appendix. 



QHAJTBR III 


G-Y CHARAOTBRISTICS OP HYPERABRUPT ..JU EG.T IONS 


3.1 Derivation of tlie C-V characteristic: 

The G-Y characteristic of the hyperabrupt junction 
can be derived from eqn, (2.7). of the last chapter which 
is again writt-en below. 


7(¥) =V„ + Y^ = 




a 


-.No.D' 


[(1 - ( 1 + ^ )e 


-¥/l 


f 1 ^ ^ -iL - 

O 0 


jr 

2D^ 


(3*1) 


This equation can be written in the following 

form. 

-l/o 

V(¥) = V [ ( 1-R.).( l-(l+l/0)e ) + ] (3.2) 

° ° 2q^ 

where, 

is equal to N^/N^ 

’C* is the normalized capacitance and is equal to 

l/¥. 

N 1^ 

’V * is equal to q, 

^ e 

A dimensional analysis of different terms of the 

eq. (3.2) reveals that *R* and O'’ are dimensionless 

c 

I I 

quantities while Y^ has the dimensions of voltage. Since 
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f t ^ 

is constant for a j unction , it can be used as a normali- 
zing factor for V(¥) . ' $he eq. (3»2) therefore, can be 
further written as, 

(1-R ) ( 1-(1+ ^ ) e (3.5) 

c 0 2 


where 



Jill. 


Y. 


is the normalized -voltage. 


3.2 Dependence of 0-7 index (21) upon impurity profile 
parameters of retrograded regions 

Ihe equation for 0-7 index can be derived from the 
0-7 relationship of hyperabrupt junction i.e. eq., (3.3). 
Differentiating both the sides of eq. (3.3)* with respect to 
* 0’ one obtains. 


d7n 
dO “ 


-l/'O 

(l,Rc)e 

0 ^ 


(3.4) 


Hence, 


n 


d (In «) 
d(ln 7k) 



[(1-R^) (l-d+l/Oe ) 4^] 
( l-R^ ) e“^/ ^ 

(3.5) 


• • ^ 

!Ehe minus sign on right hand side indicates that 

the plot of In 0 versus In 7^^ has a negative -slope. Equation 

(:>.5) also shows that n can be expressed as a function of twb 

variables, 'G' and 'R * , VJand 'n* have been calculated for 

o 

different values of ’0’ and ^R.*' using eqs. (3.3) and (3.5) 
respectively, ’n' has been plotted as a function of ’C’ 
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and '7^ in Figs. (5)and (6) respectively with*R' as 
parameter. 

It is seen from Fig, 5 that for a given value of 
‘n’ first increases reaches a maximum and then starts 
decreasing with increasing ’0’ [Fig, (5)]* For very low as 
well as for large values of "C’^'n* is seen to approach 
the value 0,5 > which is the 0-7 index of an abrupt 
junction* The dotted curve is the plot of ri ^ , as a 
function of *0*. The nature of the n versus ’0^ curves > 
of Pig, 5 can be explained as follows. 

For large values of ‘O* i.e. low values of ’W’ , the 
space charge is confined to the region near the junction# 
where impurity concentration is ver^?- high. In this 

» I 

situation# the change ^ If in W, caused by a differential 
increase Al in V#' is very low* The impurity concentra- 
tion at X =¥ and at x = ¥ + # therefore# are almost 

equal and junction behaves more or less like an abrupt 
junction , As the space charge moves further in the 
retrograded region, the impurity concentration decreases and 
same amount of change in voltage causes higher and hi^er 
change in # thereby increasing the value of ’n^ '. At 
very low values of '0* far away to the left of the 
maximum points where the space charge has covered practically 
whole of the retrograded region and has entered the 
substrate region, the junction again exhibits abrupt 
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junction characteristics, because the impurity concentration 
in substrate region does not change with the distance. 

A. careful inspection of GMg. (5) also reveals that 
maximum value of n, n^„,^ >' depends on -R ' only. A plot of 
%ax.' ® function of R^ is shown in Rig, (6). It is seen 

that ’ increases with decreasing ’'R^', This is expected 

because as R increases and approaches unity the junction 
tends to become an abrupt junction, 

C-7 index (n) has also been plotted as a function of 
’7’ [ Pig. (7)]. It is interesting to note that ’n’ is always 
maximum when is approximately unity. 

5.3 Relationships between’n ' 'R * and ’O': 

#0 

Empirical relationship expressing ’ C’ as a function of 

’n^„^ * and ’n^ as a function of R^ 'have been eobtained 
max . max c 

by linear regression method [22], Former relationship is 
given as, 

In(lOC) = 4.43 0.391n (10 - 2.8 f In (10 (3.6) 

(1.07) ( 0 . 3 ) * (1.15) 

For the- equation, R^ = 0.986 and R^ = 0.98 

The latter relationship is, as given below. 

In (10^u^„^ ) = 7.98 - 0.68 In (10 R^) + 5.29 exp (R^) (3.7) 

( 0 . 61 ) ( 0 . 012 ) (0.64) 

The values of R^ and R^ for this eq. are, R^=0,999 and 
R^ = 0.999. 

The quantities written below the coefficients, in 
brackets, are standard errors of respective co-efficients . 






# 1 ^ 


,v='> 'HH"-> ' *-^'’- 
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interpretation of and has "been given in appendix. 
Hie above expressions as well as the empirical 
expression given in Chapter II have been obtained by using 
a standard regression program available in Computer Centre 
Library, i.I.IC.Eanpur. 
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EXPB&IMBHPAL WORK 

.4.1 Fabrication Technologys 

Several techniques have been reported for fabricating 
the hyperabrupt junctions. Major among them are (i) alloy 
diffusion [l-4] (ii) double diffusion [24] (iii) epitaxial 
growth [25] and (iv) ion implantation [26]. In our work the 
double diffusion process has been used. 

4.1.1 Design Considerations: 

Since the aim of the fabrication of the hyperabrupt 
junction has been to verify the theoretically calculated results 
the device structure and the values of impurity profile 
parameters are to be decided from the point of view of verifi- 
cation of theoretical results only, 

The fabricated junction has to be a p -n type, because 
all the calculations for avalanche breakdown voltages have been 
made for p'^-n junctions only. Also for the measurement of 
true breakdown voltage of the junction, it is necessary that 
the structure of the device is mesa type, as this structure 
eliminates the effect of curvature of the junction (curvature 
lowers the breakdown voltage ) on the breakdown voltage [27]. 

Further as explained in Chapter II , the junctions with 

F 2 z 10^'^ and I 1,06 micron, behave like the abrupt 
o 

junctions. The fabrication of the junctions having' and ' L' 
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in this range, therefore, can he dispensed with* 

4*1*2 Jkhri cation li'ocess^; 

The hyperahrapt p'^-n junctions were fabricated on 
n-type silicon wafers, polished on one side* Ihe resistivity 
of the wafers was . specif ied to he in the range 3-7 ohms. cm. 
and the measured thickness was 8 mils, Hhe complete process of 
fabrication comprises of following steps. 

Step I Gleaning: 

After degreasing, the wafer was cleaned by boiling it 

in the sulnit solution ( Cone* H2S0^ + Gone* HNO^, 1: 1 by 

volume) for 10 minutes* A thorou^ rinsing in distilled 

de-ionized water (DDIW) and drying in hot air followed the 

cleaning in sulnit solution* Ihe wafer was then transferred 

to the phosphorous diffusion system* 

Step II Ihosphorous Diffusion: 

A liquid dopant POCfl^ was used as source for phosphorous 

diffusion. Beposition was done at 800® G, for 10 minutes. Ihe 

flow rates of oxygen, source nitorgen ( used to bubble the 

POOl^ ) and carrier nitrogen were maintained at 40 c.c*/min, , 

15 c.c./min. and 350 cc /min* respectively* This was followed 

by a drive -in at 1100®C for 2 hours. The calculated surface 

concentration and, diffusion length under these diffusion 

18 

conditions are 7.9x10 cm and 1,06 micron respectively* 

Step III formation of Exponentially Retrograded Region: 

Since the impurity distribution expected from the above 
described diffusion is Gaussian function , only the tail of 

• err> 25 . 
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which (about two decades below the surface concentration) can 

be closely approximated to the exponential profile, it was 

considered necessary to etch out a small thickness of wafer 

so a s to expose the surface ^ere the impurity concentration 

of original Gaussian distribution has fallen by two orders of 

magnitude . Ihis etching also gives the added advantage 

that it provides the desired low values of without 

increasing the drive in time and hence the diffusion length 

’L^' and also the characteristic length ’I', 

Ihe etchant used was a solution of HP, and 

OOOH in ratio of lt8;l by volume. Ihe etch rate of this 

solution was found to be 5.2 micron/min. illVo phosphorous 

diffused samples ( 1 and 2) were etched in this soluticn 

for 25 sec. and 30 sec. respectively. Pdr first case, the 

17 —5 

expected surface concentration was 1.03x10 cm and for the 

16 —’5 

second it was 2x10 cm 

Diffusion Length(D^)j The diffusion length (1^) is given as, 

= 2 V’D.t 

where, D = Diffusion co-efficient of the impurity atoms 
( cm, sec . ) 

t = Diffusion time (sec.) 

1 » 

This diffusion length,!^, is a parameter of Gaussian 
impurity distribution obtained from drive-in of impurities and 
is different from characteristic length, ’L’, defined earlier 
in Chapter II and III, which is related with the exponential 
distribution, an approximation to the tail of the impurity 
distribution obtained by the diffusion process. 



24 


St ep 17 Boron Diffusion: 

Boron Nitride (BN) wafers were used as a solid state 
source for boron diffusion. Before the diffttsion was carried 
out, the cleaned BN wafer was oxidised at 900° (T for 20 minutes 
.with oxygen flow rate maintained at- 200 c.c./min,- The samples 
were then loaded for p-type diffusion. Since the irnctions were 
required to be shallow .p'*~-- n type, oh*’.’’ deposition t: as carried 
out at, ' 950°C for 30 min. with nitrogen flow rate of 550 c.c./mi;.. 
The expected junction' depth ^d the mRorured sheet resistance 
were 0,26 micron and 9 ohms/ respectively . From these data 

the surface concentration . was det ''rnir.ed to be 4x10 cm. [2 3]" 

The shallow junction could also be obtained by carrying 
out deposition at high temperatures [ 1000°C-1200°C] for a 
less time duration . But' it was observed that it gives rise 
to the formation of some compound of silicon and boron t^hich 
t^.n) 0 t possible to remove. Therefore, it was decided to carry 
out the deposition at lower temperature. 

Step V Ohmic Contacts: 

The n surface i.e. the unpolished back side of the 
wafer, was lapped with 9 micron 0^ powder so as to remove- 
the p-diftu:sed layer completely-. Ohmic contact, to this side, 
then , was made by electroless nickle plating [30]. The 
front side y( P'*^- side) , being very heavily doped, is expected to 
provido Q. good ohmic coniac'b *fco iho moiai probes* 
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Step VI Mesa Etching: 

A thick solution of apiezon wax ("black wax ) was made 
in tri— chloroethylene and the small dots of this solution 
were put on the p"^- side of the wafer with the help of a 
sy^tnge* The unpolished side was completely covered with 
apiezon wax. Since the apiezon wax is not reacted upon 
the etching solution, the etching of this masked wafer 
resulted in the mesa diodes of area nearly equal to the 
area of the dots on the polished side. The wax was then 
rmmoved by dissolving it in tri-chloro ethylene. 

4.2 Measurements: 

Out of a number of junctions fabricated , a small 
number was found to exhibit a sharp breakdown. The 
measurements were taken on three representative junctions 
( in. 1 on sample 1 and jn, to and 5 on sample 2 ) 
only. The breakdown voltages were measured directly from 
the current - voltage characteristic displayed on the curve 
tracer. 

The capacitance measurements were taken on a Boonton 
Electronics Corporation, type 74C-S18 capacitance bridge 
This is a fixed frequency ( 100 KC ) bridge which reads 
an accuracy of 0,25 percent in absence of any shunt 
conductance. The respective voltages were measured on a 
’Sensitive DC Meter, Model 95A,’ manufactured by Boonton 
Electronics Corporation, The observed capacitances for the 
three junctions are plotted as a function of the voltage. 
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on log-log graph paper, in M.g. (8) a and Fig. (9)a, 

4.3 Determination of Impurity Profile in Retrograded Region: 

Impurity distribution in the retrograded region has 

been udetermined by the 0-V characteristic of the junction. 

Assuming that (i) junction is one-sided (ii) all the impurity 
near 

atoms /L the junction are completely ionized, and (iii) 
all the voltage drop occurs in the . space charge region, the 
impurity concentration at" any point can be given, in terms 
of capacitance and voltagev by the following expressions [ 31 ] »■ 

dV'. 

N(x) = — — . ■ ^ (4.1) 

/f2 dOm 

1 .q.e T 

and, X = . (4.2) 

■ ■■ 

= SB ^tal cpapacitance of the junction • 
where A = Area of the diode 

X = Distance from the junction. 

^■(x)* Impurity concentration atcith^ distance ^x’ 

•from the junction 
7 = Voltage across the junction. 

. i 

Phe eq. (4.1) can also be wri-fcten as, 

• ' V 

IT(x), = — . (4.3) 

A\q,e , n. • 

where 'n' has its usual meaning. 

fhe values of 'n’ were measured at different points of 
C-? characteristics [Fig, (8)a and Fig, (9)a)] and the values 




IMPURITY PROFILES. 
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of 'x* and ’lT(x)' were calculated by using eqs. (4.2) and 
( 4 . 3 ) respectively. ’H’(x)’ has been plotted as a function 
of 'x' (,Fig. (8)t »Fig. (9)b). fhe impurity profiles in 
retrograded regions of junction 2 and 3 are same, as 
expected. 

4.4 Results and Discussion: 

From the impurity profiles of retrograded region, 
obtained from the C-V measurements of the junction, the 
characteristic parameteirs of the impurity distribution 
for various junctions have been determined. From Fig, 

(8)b and Fig, (9)b. it is seen that* ln[F( x) ]* versus ’x’ 
plot is linear near the junction, indicating that impurity 
distribution is exponential in this region, The cross over 
impurity concentration, has been determined by extra- 

polating the strai^t line portion of the plot to x=0. 

The values of for the junctions of sample # 2 is 

observed to be lower than that for junction of sample #■ 1, 
which is expected, for the sample 2 has been etched for 
a longer time. The distance, in which (F^-Rg), which 
is nearly equal to in these junctions, drops to • l/e’ 

of its initial value, gives the value of ‘It’ ... The substrate 
impurity concentration is found to be 10 cm which 
corresponds to a substrate resistivity of 3*2 ohm. cm. This 
value is in the range of resistivity specified by the 
man^acturer (3-7 ohm cm.). 
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Both the impurity profiles, shoxm in Pig, (8)b and 
Big< have a peculiar feature and that is, the dip in the 

imphfity profile, just preceding the region of constant 
impurity concentration of substrate. Bo plausible explai- 
natioh' Caii be given at present for this unexpected result. 

The ’Via' and ’n^ * as measured for the various 

junctions^ along with the respective parameters characteriz- 
ing the impurity distribution in retrograded region are 
given in table 1« ’> given in the table, is the 

mxii p 

impurity concentration at the lowest point in the dip, 
mentioned before* 

Since the eq* (2.11), for 'Vg' does not hold 
good for the values of '1' below 0.5 micron, the value of 

calculated for these junctions. 

However, it is seen that the junctions with the lower 
have a higher breakdown voltage, which is expected (see 
Table 1). 

The comparison of theoretically calculated 

and measured values of 'n * (Table 1) brings out certain 

interesting points. In case of junctions of sample 2 , 

the measured value of’n^ , is in close agreement with the 

calculated 'n ’ [using eq. (3*7)], if in the calculation of 
max 

’n is taken to be equal to 'B . /B , But the 

max C min. o 

calculated and measured values of ’n " come out to be 
considerably different if ’Rq* is taken to be equa,l to 
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'E-ri/'B » Exactly reverse is the case for the junction . 

±5 0 

If * *' is calculated with 'Rq' eq.ual oto Eg/R^, it gives 
a value very close to the measured value of ’n _ ' while 
the agreement between calculated ' aiid measured 

is very poor if ’Rq” is given a value eq.ual to ' ^min.'^^o'* 



GONCLUSIOM 


AvalaxLclie breakdown voltage for exponentially 
retrograded silicon p-n junctions have been calculated for 
various combination of parameters of exponential impurity 
distribution, The empirical relationships [eg, (2, 11) and 
eg. v2,l2;_,, --h; :r_ volt, rf -'■ vc^ 

charge width at brealrdown as a function of the impurity 
profile parameters, have been established, The calculated i 

results mark the upper bound of breakdown voltage. In | 

actual junctions, the breakdown voltage is expected to be i 

less than this bound because of microplasmas and the ■ | 

junction curvature. 

The calculated breakdown voltages have been compared : 

with those reported by M.Shinoda [5]. It has been observed | 

that for the same impurity profiles, the breakdown voltages | 

reported by M.Shinoda are higher. However, our results for | 

breakdown voltages are expected to be in better agreement with | 

i 

the experimental results than those reported by M.Shinoda,, j 

for the reason that more accurate values of ionization rates fcr | 

i 

I 

for holes and electrons have been used for calculation of | 

I 

breakdown voltage, j 

The dependence of 0-1 index on the impurity profile 
parameter has also been studied. It is observed that nearly 

j 

in all the cases considered, the junction shows up its | 

hyper abrupt nature only in the range, 0,1cm £ C $ 10cm | 
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[Fig. (5)] where '0^' is the value of normalized capacitance 
at n = Outside this range of *0' , the junction 

approaches an abrupt junction . Moreover the peaks of the 
curves in Fig. (5) can be said to be the most important 
points because these are the points which make the hyperabrupt 
junction superior to the other types of junctions. Also it 
is only at these points that ’n' ( which is equal to n . 
at the peak) remains constant over a small range of voltages 
[ Fig. (7)]i x^hich is a requirement in moat of the applica- 
tions of varactors. The empirical relationships, relating 
*n * ’C ' and have been obtained [eqs. (3«6) and 
(3.7)] by using linear regrassion technique. 

The four equations [(2 ,11) , (2 .12 ) , (3^6) and (3.7)] 

form a part of the set of design equations for silicon 

hyperabrupt junctions. In the design of varactors, two 

different types of situations can be encountered. First is 

where the range of variation of capacitance is of primary 

interest and second xfhere the value of 'n’ is the main 

criterian of design. The starting point of the design in 

both these situations, is the assumption of an impurity 

profile .based on the requirements of cut off frequency of 

the varactor. For a given impurity profile, the above 

used 

empirical equations may be^in the design of varactors as 
explained below. 



32 


Considering the first type of situations, the minimum 
capacitance corresponds to the space charge width* wl can he 

±5 

determined by substituting and 'Ng* in equation 

(2.12). The corresponding voltage can be calculated by 
Eq. (2.11). Eor the calculation of mazimum available 
capacitance none of the four equations is helpful. It cam 
however, be determined by using Lawrence and Warner's 
plots [32], 

For the second type of situation where C-V indez 
(n) is the design criterian, 'n^ is calculated first by 

1113, X 

eq. (5*7). For this value of , ' G^* and hence 

^ maz m . m 

( = ) is then determined by eq. (3*6) . 'Wg* is now, 

m 

calculated by substituting », 'L' and 'IL,' in eq. (2,12), 

0 15 

Since in most of the applications, it is desired that the 

junction is operated at n = n , , the *W * should be 

maz m 

smaller than '¥^'. By the iteration of this calculation 

15 

procedure, a compromise between cut-off frequency, which 
is partly determined by 'E ','L' and and *n can 

O -D II13X 

be reached. 

In an attempt to verify some of the results of our 

calculations, hyperabrupt junctions were fabricated by 

double diffusion technique . The main problem encountered 

in the fabrication, was that of controlling the cross over 

concentration , 'E ' . The surface concentrations of 

o 

impurities obtained by phosphorous diffusion x^ere not in 

I 

the acceptable range of Eo*. So the surface concentration 
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was reduced iDy etching out a small thickness of wafers, ■ 

thus, exposing the surface of lower impurity concentration. 

For a given etching rate, the thus could he controlled 

by controlling the etching time only. The etchant used. 

for the purpose is felt to be having rather a high etching 

rate, resulting in a some what poor control over ’N M 

0 

attempt to decrease the etching rate of this solution by 
increasing the moderator (CH^COOH) was made. 

No marked decrease in etching rate was observed, even when 
the composition of etchant was changed from 1 part HP + 8 
parts HNO^ + 1 part CH^OOOH to 1 part HP + 8 parts- HNO^+ 3 
parts CH^COOH. Further increase in CH^COOH content was 
found to give rise to a thin whitish hazy deposition on the 
surface. Por the improvement in control of 'N^' , it is 
essential that an etchant of much lower etching rate, than 
that used in the fabrication, is developed. Por mesa 
etching, the photolithography should certainly have given 
better results as it gives the diodes of well defined 
geometry , thus making it possible to measure the area of 
the junction accurately. The evaporated gold dots probably 
can also be used as a replacement to blackwax dots, which 
were used for the mesa etching in the fabrication of hyper- 
abrupt junctions. 

The theoritical results could not be verified 
experimentally. The experimental results, obtained farom 
the measurements on a few hyper abrupt junctions fabricated. 
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are rather inconsistent with the theoretical results^ 
However, heoause of the limited number of observations it is 
not possible to draw any definite conclusions. 


/ 



APPENDIX 


The method of regression analysis consists of postula- 
ting a relationship between a dependent variable an,d independent 
variables and then determining the co-efficients of the 
independent variables, for the least sum of squares of residuals 
[ residual is defined as the defference in the observed and 

estimated value of dependent variable]. The quantities -R^ 'and 
“ 2 * 

R help in assessing the goodness of the equation in estimating 
the value of dependent variable for a given set of values of 
independent variables. 

(i) R^ j R^ is given by 

R^ = 1 — 

(Y-Y)" 

where Y = observed values of dependent variable 
Y = mean value of dependent variable 
e^ = sum of square of residuals. 

It is obvious from the expression for R^ that the value 
of *R^ 'approaches unity when the value of residuals, *e* 
approaches zero, i.e. when the value; of dependent variable 
estimated by the equation, approaches the corresponding 
observed value of dependent variable for each observation , 

(ii) R^ : This is given by the relation, 

V(e) 

R^ = 1 - — 

V(Y) 



where 


■^(e) = Variance of residual = 

ii rt 

' f 

v;(y) = Variance of dependent variable 
= (Y-Y)^/(n-l) 

= Degrees of freedom 

K ■ = Total number of co-efficients including 

constant term 

n ' = Total number of observations. 

With, the addition of variables '/terms to the equation/ 
necessarily increases regardless hf the relevance or 
otherwise of the indluded term. This ia beipause of decline 
in • e^‘. Here in, ,the value of assumes importance. 

Dor the degrees, of freedom also decreases and hence 'V(e)’ 
may become less or more, depending upon the relative 
r edubtion in ’ e^ * and • Thus / if the ad.dition 

' of a' increases also, it implies that after the 
inclusion of term, . the expression has a better power of 
prediction, due to the lessening of the variance of the error 
of prediction. 
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